Cosmological CPT violating effect on CMB polarization 
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A dark energy scalar (or a function of the Ricci scalar) coupled with the derivative to the matter 
fields will violate the CPT symmetry during the expansion of the Universe. This type of cosmological 
CPT violation helps to generate the baryon number asymmetry and gives rise to the rotation of 
the photon polarization which can be measured in the astrophysical and cosmological observations, 
especially the experiments of the cosmic microwave background radiation. In this paper, we derive 
the rotation angle in a fully general relativistic way and present the rotation formulas used for the 
cosmic microwave background data analysis. Our formulas include the corrections from the spatial 
fluctuations of the scalar field. We also estimate the magnitude of these corrections in a class of 
dynamical dark energy models for quintessential baryo/leptogenesis. 
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I. INTRODUCTION 



In the standard model of particle physics, CPT is a fundamental symmetry. Probing its violation is an important 
way to search for new physics beyond the standard model. Up to now, CPT symmetry has passed a number of high 
precision experimental tests in the ground-based laboratory and no definite signal of its violation has been observed 
[l[. So, the present CPT violating effects, if they exist, should be very small to be amenable to the experimental 
limits. 

The CPT symmetry could be violated dynamically in the expanding Universe. To show it, consider a scalar boson 
<f> which effectively couples to a fermion current J M , with the Lagrangian given by 

Ant = ^V^J" , (1) 

where c is a dimensionless constant and M is the cut-off scale. The interaction in is CPT conserved; however, 
during the expansion of the Universe, the background value <f> does not vanish and CPT is broken spontaneously. 
This type of CPT violation occurs naturally in theories of dynamical dark energy and has interesting implications in 
particle physics and cosmology. In models of quintessential baryo/leptogenesis [3, 3|, kyl , the scalar field <f> in JT]) is the 
dark energy scalar (quintessence [ESQ, fc-essence phantom [9|, quintom [xol. Till] etc.). In the early Universe the 
field (f> with the interaction in ([1]) generates the baryon number asymmetry and at late times it drives the accelerating 
expansion of the Universe. One of the features of these models is a unified description of the present accelerating 
expansion and the generation of the matter and antimatter asymmetry of our Universe. Furthermore, differing from 
the original proposal for spontaneous baryogenesis by Cohen and Kaplan [T^], since the dark energy scalar has been 
existing up to the present epoch, the corresponding CPT violation could be tested in laboratory experiments and 
cosmology. In Refs. 0, H|, we have pointed out that, to produce the enough baryon number asymmetry, the dark 
energy should be significant in the radiation-dominated epoch. This is the case if the dark energy has the tracking 
behavior, i.e., its density decays almost at the same rate with that of radiation as the Universe expanding. Along 
this line, the gravitational baryo/leptogenesis [H, [l4[ has been proposed in which a function of curvature scalar R 
replaces the <f> field in JT]). There are other motivations in the literature, e.g., Refs. 

mmmmmnii 

The current J M in Eq. (Q]) is not necessary to be the baryon current for baryogenesis. It could be other currents 
which are not orthogonal to Jg or J I E \_ L - In [21| , we have proposed, for example, that J M is the left-handed part of 
the B — L current J^ B _ L ^ L ■ Besides the generation of baryon number asymmetry, this kind of coupling will bring a 
new effect to the photon sector. This is because Jm—L) L 1S anoma l° us under the electromagnetic interaction 

V V J (B-L) L 37r r V* r ■ V 1 ) 
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Hence the interaction in Eq. (JTJ) would induce the following effective coupling through the anomaly equation: 

Ant = ™f^f = V»A V F» V , (3) 

where A v is the electromagnetic vector potential, F M „ = V ^Ay — V ' is the strength tensor and F 1 *" — l/2e tlvf " T F pa 
is its dual. This Chcrn-Simons term leads to the rotations of the polarization vectors of photons when propagating 
over the cosmological distance [22| 1 . The change in the position angle of the polarization plane Ax, characterizing 
CPT violation in this scenario, can be obtained by observing polarized radiation from distant sources such as radio 
galaxies, quasars [13, and the cosmic microwave background (CMB) [28|, [29J. Assuming the rotation angle is 
homogeneous and isotropic Ax — Ax, the CMB power spectra would be rotated as [29l. [30j 
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In the formulas above, the quantities with the superscript obs are those observed after the rotation. T, E and B 
represent the temperature, the electriclike and magneticlike polarization modes, respectively. 

In Ref. j29j, with Feng and Li, we did the simulations on the measurement of Ax with the future high precision 
CMB experiments, CMBPol [3l[ and PLANCK j$2| using the rotation formulas We pointed out that in such 
experiments the EB spectrum will be the most sensitive probe of such CPT violation, this is because the EB power 
spectrum is generated by the rotation of the EE power spectrum, which is a more sensitive probe of the primordial 
fluctuations than the TT and TE spectra. In [301 ] . with Feng, Xia, and Chen, we first found that a nonzero rotation 
angle Ax = —6.0 ± 4.0 deg (lcr) is mildly favored by the CMB polarization data from the three- year Wilkinson 
Microwave Anisotropy Probe (WMAP3) observations [HI, 
of BOOMERanG (hereafter B03)[H, M, IU (see also Ref 
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|36l, 13711 and the January 2003 Antarctic flight 
43 . l44j). This is a signal in some sense of the 
cosmological CPT violation mentioned above. Later on, Cabella, Natoli, and Silk [45[ performed a wavelet analysis 
of the temperature and polarization maps of the CMB delivered by WMAP3. They set a limit on the rotation angle 
A% = —2.5 ± 3.0 deg (lcr). This is consistent with our result because they considered WMAP3 data only. Using the 
full data of B03 and the WMAP3 angular power spectra, one of the authors (X.Z.) with Xia et al. 14611 has found 
that Ax = —6.2 ± 3.8 deg (lcr). This result improved the measurement given by our previous paper [301 ] . Recently, 
the WMAP experiment has published the five-year results for the CMB angular power spectra which include the TB 
and EB information [47|,[48[. They used the polarization power spectra of WMAP5, TE/TB (2 < I < 450) and 
EE/BB/EB (2 < I < 23) to determine this rotation angle [49| and found that Ax = -1.7 ± 2.1 deg (lcr). However, 
when B03 data are included, one of the authors (X.Z.) with Xia et al. found that Ax = —2.6 ±1.9 deg (lcr). 
Again, a small CPT-violating effect is mildly detected by current data. 

We note that the rotation formulas given in ^ are valid only for a homogeneous and isotropic rotation angle 
and are obtained in the Minkowski spacctimc or the spatially flat Friedmann-Robertson- Walker spacetime which is 
conformally equivalent to the former. This is expected to be a good approximation when the coupled scalar field (j) is 
the dark energy or the function of the curvature scalar because in these cases 4> is very homogeneous in the observed 
Universe, while the accompanied perturbations are much smaller. Usually, its background part makes the dominant 
contributions. One of the aims of this paper is to study the secondary effects due to the perturbations of (j), which leads 
to the anisotropics of the rotation angle. For this purpose, we first study the Maxwell theory modified by the Chern- 
Simons term in the general curved spacetime and investigate the possibility of obtaining the rotation angle in a fully 
general relativistic framework. The spatial fluctuations of the scalar field make the rotation angle inhomogeneous and 
anisotropic and bring higher order corrections to the rotation formulas Eq. (01 . Specifically, we evaluate the magnitude 
of these corrections in the models of tracking dark energy, as required by the quintessential baryo/leptogenesis, and 
found the corrections are negligible. However, for some other models, the corrections could be sizable. This paper is 



1 Such a term also breaks the Einstein equivalence principle l23t 1 24tl in the short wavelength limit and breaks causality in the long 
wavelength limit [25l, |26H . 
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organized as follows. In Sec. II, we present the relevant equations of the modified electromagnetic theory under the 
geometric optics approximation. In Sec. Ill, we study the generalized Stokes parameters and the changes in CMB 
power spectra in Sec. IV. In Sec. V, we evaluate the corrections due to quintessence fluctuations in the quintessential 
baryogenesis model. Section VI is the conclusion. 

II. BASIC EQUATIONS 

The full Lagrangian of the Maxwell theory modified by the Chern-Simons term (J3J) (without other sources) is 

C = -~F flv F^+p fl A v F>" . (5) 

This Lagrangian is not gauge-invariant, but the action integral S — J Cd 4 x is gauge-independent because p p is defined 
in (|3|) as the derivative of the scalar field. The equation of motion can be obtained through varying this Lagrangian 
with respect to A v : 

V M F"" = -2 Vil W v . (6) 

The right-hand side of the above equation is brought by the Chern-Simons term. But the identity is unchanged: 

V p F pCT + V p F CTp + V CT F pp = . (7) 

We will study these equations in a gauge-independent way though it is easier to do it by choosing the Lorenz gauge 
(2lj . For this purpose we make a differentiation to Eq. ||7J) and get 

□F pCT + 2V p {p p F^) - 2V a (p p F^ p ) - [F a p R aa - F a a R ap - F^ a R appa } = , (8) 

where R aa and R alipa are Ricci and Riemann tensors, respectively. 

Since we are studying the light which propagates in the cosmological scales, the geometric optics approximation 
(GOA) applies very well. With this approximation, the solution to the equation of motion is supposed to be 

where we made a complexification to the electromagnetic field, but e is a small real parameter and S is a real function. 
This ansatz means that the phase of the wave varies much faster than the amplitude. We define the wave vector as 

K = v p s , (io) 

which represents the travel direction of the photon. 

Substituting the ansatz ((9j into Eqs. ([8|) and ([7]) and dropping out the terms containing Ricci and Riemann tensors, 
we have 

2i i 1 

□ («V + fbpo + .-.) + — fc M V AI (a pCT + eb pa + ...) + -{V p k^){a pa + eb pa + ...) - -^k^{a pa + eh pa + ...) 

= -2[(V pP^ia^ + eb pa + ...) + ^(V p a pCT + eV p 6 M<T + ...) + ^^(o^ + eb pa + ...)] + 2[p «-» a] (11) 

and 

i 

[V p (a pCT + tb pa + ...) + -k p (a pa . + eb pa + ...)] + [pafi] + [afip] = . (12) 
At the leading order of the GOA, Eq. (THJ) gives 

(13) 

which implies that a pa should have the following antisymmetric form: 

a P a = k p a a - k a a p . (14) 
Then we collect the terms of Eq. ([i"Tj) at the orders of 1/e 2 and 1/e, respectively. At the order of 1/e 2 , we have 

k^ = . (15) 
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The propagation equation of W can be obtained by differentiating the above equation again: 

= S/ V (k^) = 2V M SV l/ V M S' = 2V Al SV M V„S' = 2fc M V A1 fc 1/ . (16) 

This is a geodesic equation. The vector k^ defines an affine parameter A which measures the distance along the light 
ray: 

= < 17 > 

We can see from (ITBl that is parallelly transported along the light curve x p (A). In other words, photons travel 
along null geodesies. These results are the same as those of the standard Maxwell theory. The modification due to 
the Chern-Simons term appears at the order of 1/e: 

Va" + °-a v = -V^ v <>°k p a a , (18) 

where we have considered Eq. (TT4"]) and defined the operator T> = fc M V M . The quantity 9 — V^fc^ describes the 
expansion of the bundle of the light. Without the modification, the right-hand side of the above equation would 
vanish. Its physical meaning is that the polarization vector of the photon is not parallelly transported along the light 
ray as we will see in the next section. In addition, by applying the GOA to the original equation 

V„F"" = -2p^ v (19) 

we have 

k^ = . (20) 
The basic results we got above are Eqs. (fT5)l and (TIB)) with two orthogonality relations (fl"5|) and (f^O]) . 



III. STOKES PARAMETERS 



It is convenient to use the Stokes parameters to study the polarization of radiation. The four Stokes parameters 
are well defined in Minkowski spacetime (the inertial frame). Considering a monochromatic electromagnetic wave of 
frequency ujq propagating in the +z direction 

E x = a x (t) exp[i(uj t - d x (t))], E y = a y {t) exp[i(uj Q t - 6 y (t))}, (21) 

the Stokes parameters are defined as the time averages 

/ = (e x e* x ) + (e v e;) , 
Q ee (E X E*) - (e v e;) , 
U (E X E;) + (E*E y ) , 

V ee i[(E x E;) - {E*E y )}. (22) 

In general relativity, these definitions should be generalized. This can be done by using the tetrad formalism. 
A tetrad is a set of four orthogonal unit basis vectors e r a y with a = 0,1,2,3. At each point x, we can attach a 
tetrad which transforms between the coordinate frame and the local inertial frame at x. For a vector field B^(x), its 
components in the local inertial frame are 

Ba = e^B^ . (23) 

The Latin indices are lowered and raised by the Minkowski metric r] ab , the Greek indices, however, by the coordinate 
metric g^ v . The tetrad has the following properties: 

9^\ a f\ b ) = Vab , V ab e\ af \ b) = <T . (24) 

We can set the tetrad frame at each point as follows. Consider the rest frame of the free fall observer, in which the 
four- velocity is u a = Sq. Furthermore, we require the observer to see the light traveling along the +z direction, and 
hence k a = uj(Sq + 5f). So, after transforming to the coordinate frame, 



U — L (a) " — e (0) ' 



(25) 
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and 

W = e\ a) k a =u>{u» + e» {3) ) . (26) 

Hence, 



e" (0) =u", e^ 3) = i(^ (27) 



where u> = k^u^ 1 is the frequency measured by the observer. The other tetrad vectors and e^j are unit spacelike, 
orthogonal to each other and to e ^ y ^(3)! an d therefore orthogonal to fc M . 

The electric vector in general spacetime for an observer with four-velocity is defined as 

e p. = F nu Uv _ (23^ 

At the leading order of the GOA mentioned at the last section, it is 

£7" = a^^e^' = (Jfe"a" - k v a^)u v e iS/t . (29) 
Transforming it to the local inertial frame, we get the x and y components of the electric field in this frame easily: 

E X = E 1 = E^ (1) , E y = E 2 = E^ {2) . (30) 

In the local inertial frame, the definitions of the Stokes parameters (|2"2"|) are applicable. By applying the above 
equations to (|2"2"j) , we get the general expressions of the Stokes parameters in curved spacetime [Sj, |52| : 



I 



Q = c^L^e^e 1 ^ -e^ {2) e v {2) ) 
U = w%^(e^ (1) e* (2) +e ,i (2) e i/ (1) ) 

V = iu 2 L tlv {e li {1) e v (:2) - e M (2) e" (1) ) , (31) 

where L Mt/ =< a^a* > satisfies the following equation making use of Eq. (|18[) : 

2?i^ + 9L^ = - Pa k p (e al3 \L^ + e a P\L m ) . (32) 

We can see that the Stokes parameters are coordinate scalars but not Lorentz scalars. We require the tetrad frames 
to be not physically rotating. In order to do that, we set the tetrad vectors at each point so that e ^ and e M j 2 j 
are parallelly transported along the light curve. So it is straightforward to get the propagation equations of the four 
parameters along the light curve: 

VF + 0F Q = , (33) 
VF X + 6F l = 2p^F 2 , (34) 
VF 2 + 9F 2 = ~2p^F 1 , (35) 
VF Z + 6F 3 = , (36) 

where F a = (I, Q, U, V)/oj 2 . Equation (33) means the conservation of the light flux. Equation (36) indicates 
that the Stokes V, which describes the net circular polarization, vanishes if it is zero at the beginning. This is the 
case for CMB where the polarization is produced at the last scattering. Since the Stokes V cannot be produced 
by Thomson scattering, it remains zero afterwards. In short, the net circular polarization remains vanishing even 
in the presence of the Chern-Simons term. The terms in the right-hand sides of Eqs. (34) and (35) are the effects 
of the Chern-Simons term which rotates the polarization angle of the light. The polarization angle defined by 
X = 1/2 arctan (U/Q) = 1/2 arctan (£2/^1) satisfies 

V X +P^ = 0. (37) 

This angle when measured at the point / is rotated by 



Ax = Xf - Xi = - I P^dX = - I p^dx 11 

J i J i 



(A) , (38) 



6 



compared with that at the point i when the photon was emitted. From ((3|), = — (2ca em )/(37rA/)<9 M 0, the rotation 
angle is given by 



Defining 



F± = F 1 ± zF 2 , (40) 



it satisfies from Eqs. ([34]) and (|35|) that 

F£ = exp (- y exp {±i2A X ) ■ (41) 

The Chern-Simons term modifies the result by merely adding the rotation factor exp(±i2Ax). Hence the observed 
Stokes parameters should be 

(Q ± iU) obs = exp (±i2Ax)(Q ± iU) ■ (42) 

This is the basic result obtained in this section. It describes the rotation of the polarization of a single bundle of light. 
It is the starting point to study the rotated CMB power spectra in the next section. 

IV. CMB POWER SPECTRA 

In order to analyze the CMB map, we usually make multipole expansion. In the flat Universe 2 , we can expand the 
temperature and polarization anisotropics in terms of appropriate spin- weighted harmonic functions on the sky [54j : 



T(fi) = y^aT,imYim(n) 

I in 

(Q±iU){n) = J2 a ±2M±2 Y im(n) . (43) 



Irn 

The expressions for the expansion coefficients are 



a T ,i m = J dnY^ m {n)T{n) 
a±2,im = J dn ±2 Y l * m {h){Q±iU){h) . (44) 
Instead of a,2,i m and a_2.z m , it is convenient to introduce their linear combinations 

a£,Im = —(a*2,lm + a-2,2rn)/2 

O.B,lm = K a 2,lm ~ CL-2,lm)/ 2 - ( 45 ) 

The power spectra are defined as 

( a *X>Mm> a X,lm) = Cf X &V^ m 'm (46) 

with the assumption of statistical isotropy. In the equation above, X' and X denote the temperature T and the 
E and B modes of the polarization field, respectively. For Gaussian theories, the statistical properties of the CMB 
temperature/polarization map are specified fully by these six spectra. In the standard case, Cj B = Cf B = 0. 
Considering the rotation in Eq. (|4"2"]) , the expressions for the expansion coefficients become 

afhm = Jdn ± 2 Y l * m (n)(Q ± iU) obs (h) = j ' dSl ±2^(n) exp (±*2A X )(g ± iU){h) , (47) 



2 For the treatment of CMB anisotropies in open and close Universe, please see 
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and aT,im remains unchanged. The rotation angle in (|39|) depends on time as well as space generally. It can be 
separated as the background part, which is homogeneous and isotropic, and the perturbation, which is randomly 
distributed on the sky: 

A x = Ax + A6 X , (48) 

where 

9fYV 

Ax=^^[^o)-^ec)] (49) 

^X= ~^^S(f)(xdec, Vdec) ■ (50) 

In the above equations the subscript indicates the present values and dec means the values at the time of matter- 
radiation decoupling. The homogeneous part Ax is the same one that appeared in the previous rotation formulas Q. 
The final value of the fluctuation 5(f>(xQ, rjo) is neglected because it only gives rise to a dipole contribution due to our 
motion with respect to the CMB frame. In the flat Universe, Xd ec = (vo ~ Vdec)™ when putting the observer at the 
origin of the coordinate system. Similar to the studies on Faraday rotation of the CMB polarization by a stochastic 
magnetic field (55j . we expand ASx on the sky: 

A5x = ^2h m Yi m {n) , (51) 

I in 

and define its angular power spectrum as 

(tf, m ,b lm ) = C?5 vl 8 m , m , (52) 

where we have also assumed statistical isotropy of b\ m . This angular power spectrum is related to the power spectrum 
of 5<p at time r\dec-, which can be seen from the following discussions. Expanding 5<j)(xdeci Vdec) m terms of Fourier 
functions, we have 



(x dec , Vdec) = J j^j^ 0fc(% ec )e lfc '' A ; 

d 3 k v 



(2tt) 3 / 2 

(Vdec) J2 ( kAf l)YL ih)Y lm [ft) , (53) 



i 

d 3 k 
(2tt) 3 / 2 



where A77 = 770 — rjdec, ji is the spherical Bessel function and Pi is the Legendre polynomial. Comparing it with Eqs. 
(|50)) and (|5l|), we get 

h m = ~^ri l J j0^^ k (mec)ji(kAr,)Y; m (k) . (54) 

With the help of the definition of the power spectrum of 5(f>, 

2?r 2 

(0* k '(Vdec)(t>k('ndec)) = -j^-V^k, Vdec)8 3 {k - fe') , (55) 

we can find that Cf in Eq. ([52")) is 



Jl 9ttM 
and 



C} = ^Sfc J fv^k, Vd ec)jf(kA v ) , (56) 



9irM 2 
i 

where 8<j) — Sc/)(xdec, Vdec) and the arguments Xd ec , Vdec are suppressed in the following 
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With these formulas, we can calculate the coefficients after the rotation: 
<•■!,.„ ■ ■ I dn ±2 Y* n (n)(Q±iU)° bs (n) 

= exp (±i2Ax) ^2 a ±2,hrm J dQ ±2Y* m (n) exp (±2iASx) ± 2 Yi imi (n) 



iimi 



= exp(± l 2Ax) E a ±^ mi Ft hmi • (58) 
In the last equality, we have denned 

Ft hmi = JdQ ±2Yf m (h) exp (±2*A<5 X ) ± 2 Y limi (n) . (59) 

So 

a °E%n = o E [( e ' 2A ^^ 7 imi 1 m 1 + e '^^^Imlimi ) ttE A m l + 'l 6 ' 2 ^^!^!!!!! ~ e l2A *F lmlimi )aB,l 1 m 1 ] , 
iiroi 

<i m = jE[H)( el2A ^U mi - e ~ 42A *^^ ■ ( 6 °) 



2 

iimi 



To calculate the observed correlations of T, E and J3, we make the following assumptions: (i) the rotation field 
X or <j> is uncorrelated with the primordial T, E, and -B modes; (ii) the rotation angle is small everywhere. In 
addition, we have Cj B = Cf B = for primordial modes. Hence, we need only to calculate the following six corre- 
lations: (F lml , m ,), Ylhmi C{[ X (^m'iimi-^mlimi)' Sz imi Qf^ (-^i'm'ii ^imil "H ^ ' Edroi Qf^ (^i'm'limi-^mlirai) 

and E/ imi C^ x \F^* n , lim F{ mlimi ). Up to the quadratic order of AS X , we have 

{FLi>m>) - (1 ± 2iA ^ - 2A£ x 2 )<^<5w = (1 - 2(A<S x 2 )) ( 5 i ^w , (61) 

and 

/i \ Vm'l\m\ lml\m\ ' 

- E C h X ' J dn ' dn ^ ~ 4 < A ^ 2 ) + H^ X (n')AS X (n))} 2 3W(»') 2 ^ mi (n') 2 Y^(n) 2 r, imi (n) 
= (1- 4(ASx 2 ))C xx 'S w S mm , +4 £ Qf X 'Cf 2 1 dfl'dfi a Y,, TO ,(n') 2 y^ mi (n')^^') a*&(n) 2 Y, imi (n)Y, 



hm 2 



l 1 m 1 l 2 m2 



The remaining integrals in the above equation may be expressed in terms of the Wigner-3j symbol through the general 
relation 1561] : 



So 

E ^^mx-FiUmx) - {l-4(±Sx>))Cr'Sw6 mm . +E ^ + + ^ ( 2 -2 V*"*"" ' 

hmi hh 

where we have used the orthogonality relation of the 3j symbol 

£(2l + l)( 1 h h )( V > h h )=8u>6 mm >. (65) 
^ v ' \ m mi m 2 J \ m mi m 2 / 
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Similarly, we can find that 



limi limi 



and 



EQXX' (p~* T? + \ _ (jXX'/T?+* Tp \ 

Zl \ l'm'l\ra\ lvnl\m\l / j l\ \ Vm'l\m\ Imlimi' 
limi limi 

* (1 - 4(Ag x 2 ))Cf*V W + £(-1) w Cff (Ml + 1){2h + % l *) 2 5u>5 mm ,, (67) 

where L = I + li + I2 and we have used the permutation property of the 3j symbol 

1 h h \ = {-l) L ( 1 h h ) . (68) 
-m —mi — ui2 I V m mi m,2 y 

Consequently, we obtained the rotation formulas of the power spectra 

n TT,obs _ n TT 
°Z £ ' 

c T B)0bs = C p cos( 2Ax)(l-2<A5 X 2 )), 

c ,t B ,o6 s = c^sin(2Ax)(l-2(A5 X 2 )), 

Cf B > o6s = [Crcos 2 (2Ax) + Crsin 2 (2Ax)](l-4(A ( 5x 2 )) 

1 ?! ; 2 \ 2 (2fx + l)(2 f a + l) 

2 -2 ) 2ir 



+ £( 2 -2 ) — 9; + ( ~ 1)L+1 C ° S ( 4A ^)]^ + t 1 + cos (4A*)]Cff B } , 

Cf B '° b = [Cf £ sin 2 (2Ax) + Cf s cos 2 (2Ax)](l-4(A ( 5 x 2 )) 

+ Efs -2 ) 2 ^ + ^ + ^ + c » s (4Ax)]^ g + [1 + (-1) L+1 cos (4A X )]C7**} , 



ilia 



Cf s ' obs = isin(4Ax)(Cf £ -Cf s )(l-4(A ( 5 x 2 )) 

(^Efo -2 V + of 2 + a <%(-l) J+1 (Cg* - O ■ (69) 



sm 



In comparisons with those in Eq. (4), Eq. (69) included the corrections from spatial fluctuations. 

From Eq. (69), we can see firstly that C l '° s and C l : ° s are proportional to sin (Ax), which vanish when A\ — 0. 
This is understandable because CPT is violated only by the background field. Second, we find that 

£(2/ + l)(C EE ' obs + cf B ' obs ) 



$> + l)(Cr + Cn(l-4<A*x 2 » + £( 2 \ o 2 ) (M + 1)(Ml + 1)(2b + 1) Cg(Cg J ' + 

Z ZZjia ^ ' 

£(2; + l)(Cf £ + Cf s )(l - 4(A5 X 2 )) + 4<A5x 2 ) £(2*i + 1){C EE + C?*) 
1 h 

£(2; + l)(Cf £ + Q BB ) , (70) 



I 

where we have used another orthogonality relation of the 3j symbol 

2 



V(2Z + 1) 1 11 12 =1. (71) 

' y —mi — m-i m\ rri2 J 

The equality in Eq. ([70)1 is the direct consequence of invariance of Q 2 + U 2 under the rotation 
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V. THE EVALUATION ON THE MAGNITUDE OF CORRECTIONS IN THE TRACKING DARK 

ENERGY MODEL 

Equations (|69|) indicated that the most important corrections appear at the order of (A5x 2 )- In this section, we con- 
sider a model for quantitative estimation on the corrections. Specifically, we take the quintessential baryo/leptogenesis 
model as we mentioned in the introduction. For such a model, we have 

<^X 2 > = SPW 2 > - ^W 2 ) ■ (72) 

As was pointed out in 2], to generate enough baryon number asymmetry, the quintessence field 6 should have 
tracking behavior, which happens, for example, in the Albrecht and Skordis model (57[- In the following, we will 
evaluate (A^x 2 ) in this model. We consider the perturbed metric in the Newtonian gauge: 

ds 2 = (1 + 2<P)dt 2 - a 2 (l - 2<$>)dx i dx i , (73) 

where t = J adrj is the cosmic time and $ is the gravitational potential. The linear perturbation equation of the 
quintessence is 

... V 2 

56 + 3HS6 - — 66 + V" (6)66 = 40$ - 2V'<f> . (74) 
a z 

In the above equation, the dot denotes the derivative with respect to t. The general solution to this equation is 
decomposed into two parts: the adiabatic mode and the isocurvature one. For the model of quintessence with a 
tracking solution, the isocurvature perturbation decays away quickly 0, [H| . We need only calculate the adiabatic 
perturbation, which satisfies the adiabatic condition 

% = *4, (75) 
P P 

i.e., 

-!<?>■ <™> 

With the equation above we can find that at the time of matter-radiation decoupling (in the matter-dominated epoch) 
the adiabatic perturbation of quintessence on large scales is 

66=^ = ^=^M pl ^ , (77) 

where we have considered the exact tracking behavior of quintessence, ws = w m = and the well-known result 
$ = constant in the matter-dominated era. The parameter f2^ < 10~ 2 59, [6^] is the density of quintessence at this 
time, and M p i = l/^&irG ~ 10 18 GeV is the reduced Planck mass. So, with (<E> 2 ) — 10~ 10 , we have 

M 2 , M 2 , 
<A<W 2 > - 10- 7 ^<* 2 > - 1(T 17 ^ • (78) 

If (A^x 2 ) <C 1, the corrections in Eqs. (f69|) can be neglected safely, which happens for the cut-off scale M 3> 
10~ 8 Af p ; ~ 10 10 GeV. In the quintessential baryo/leptogenesis model, M ~ 10 8 Td where To is the decoupling 
temperature of lepton number violating interaction and is around 10 11 GeV 14 1. 



VI. CONCLUSION 



In this paper, we have studied the effects of the interaction with the derivative coupling of the scalar field to photons 
given by the Chcrn-Simons term in the general curved spacetime. Under the geometric optics approximation, we have 
obtained the general form of the rotation angle in a gauge-invariant method. We have calculated the corrections 
brought by the spatial fluctuations of the scalar field to the rotation formulas. These corrections exist due to the 
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dynamics of the scalar field 3 ; however, they have not been considered in the literature on the CMB data analysis. We 
have estimated the magnitude of the corrections in a model of scalar field for the quintessential baryo/leptogcnesis 
scenario and fortunately found that the corrections are very small and can be neglected safely in the fit to the CMB 
data. The same techniques can be applied to the case of gravitational leptogenesis in which the coupled scalar is the 
function of the gravitational field. Similar techniques can be developed to other cases in which the Chern-Simons 
parameter has other origins. For example, in a more complicated case, where the parameter is not statistically 
isotropic or even has no power spectrum, the space components of in Eq. §3§ will bring the correlations between 
a>T,i'm' an d a E.im of different I'm! and Im and so on. These complications are beyond the scope of this paper, and we 
leave them in the future work. 
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